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Abstract 

Titanium dioxide is one of the most widely investigated oxides. This is due to its broad range of ap- 
plications, from catalysis to photocatalysis to photovoltaics. Despite this large interest, many of its bulk 
properties have been sparsely investigated using either experimental techniques or ab initio theory. Further, 
some of TiCVs most important properties, such as its electronic band gap, the localized character of ex- 
citons, and the localized nature of states induced by oxygen vacancies, are still under debate. We present 
a unified description of the properties of rutile and anatase phases, obtained from ab initio state of the art 
methods, ranging from density functional theory (DFT) to many body perturbation theory (MBPT) derived 
techniques. In so doing, we show how advanced computational techniques can be used to quantitatively 
describe the structural, electronic, and optical properties of Ti02 nanostructures, an area of fundamental 
importance in applied research. Indeed, we address one of the main challenges to Ti02-photocatalysis, 
namely band gap narrowing, by showing how to combine nanostructural changes with doping. With this 
aim we compare TiCVs electronic properties for OD clusters, ID nanorods, 2D layers, and 3D bulks using 
different approximations within DFT and MBPT calculations. While quantum confinement effects lead to a 
widening of the energy gap, it has been shown that substitutional doping with boron or nitrogen gives rise to 
(meta-)stable structures and the introduction of dopant and mid-gap states which effectively reduce the band 
gap. Finally, we report how ab initio methods can be applied to understand the important role of TiC>2 as 
electron-acceptor in dye-sensitized solar cells. This task is made more difficult by the hybrid organic-oxide 
structure of the involved systems. 

PACS numbers: 
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I. INTRODUCTION 



Ti0 2 has been one of the most studied oxides over the past few years. This is due to the 
broad range of applications it offers in strategic fields of scientific, technological, environmental, 
and commercial relevance. In particular, Ti0 2 surfaces and nanocrystals provide a rich variety of 
suitably tunable properties from structure to opto-electronics. Special attention has been paid to 
Ti0 2 's optical properties. This is because Ti0 2 is regarded as one of the best candidate materials 
to efficiently produce hydrogen via photocatalysis 12 . Ti0 2 nanostructures are also widely used in 
dye-sensitized solar cells, one of the most promising applications in the field of hybrid solar cells^. 

Since the first experimental formation of hydrogen by photocatalysis in the early 1980s;? Ti0 2 
has been the catalyst of choice. Reasons for this include the position of Ti0 2 's conduction band 
above the energy of hydrogen formation, the relatively long lifetime of excited electrons which 
allows them to reach the surface from the bulk, Ti0 2 's high corrosion resistance compared to other 
metal oxides, and its relatively low cosfr^. However, the large optical band gap of bulk Ti0 2 (~ 
3 eV) means that only high energy UV light may excite its electrons. This effectively blocks most 
of the photons which pierce the atmosphere, typically in the visible range, from participating in 
any bulk Ti0 2 based photocatalytic reaction. On the other hand, the difference in energy between 
excited electrons and holes, i.e. the band gap, must be large enough (> 1 .23 eV) to dissociate water 
into hydrogen and oxygen. For these reasons it is of great interest to adjust the band gap e gap of 
Ti0 2 into the range 1.23 < e gap < 2.5 eV, while maintaining the useful properties mentioned 
above 5 . 

With this aim, much research has been done on the influence of nanostructure^— and dopants 
^LIri2 on Ti0 2 photocatalytic activity. For low dimensional nanostructured materials, electrons 
and holes have to travel shorter distances to reach the surface, allowing for a shorter quasi-particle 
lifetime. However, due to quantum confinement effects, lower dimensional Ti0 2 nanostructures 
tend to have larger band gaps^. On the other hand, although doping may introduce mid-gap 
states, recent experimental studies have shown that boron and nitrogen doping of bulk Ti0 2 yields 
band gaps smaller than the threshold for water splittingi^. This suggests that low dimensional 
structures with band gaps larger than about 3.0 eV may be a better starting point for doping. 

Recently, several promising new candidate structures have been proposed 11 . These small (R < 
5 A) Ti0 2 nanotubes, with a hexagonal ABC Pt0 2 structure (HexABC), were found to be surpris- 
ingly stable, even in the boron and nitrogen doped forms. This stability may be attributed to their 
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structural similarity to bulk rutile Ti02, with the smallest nanotube having the same structure as a 
rutile nanorod. 

A further difficulty for any photocatalytic system is controlling how electrons and holes travel 
through the system — . For this reason, methods for reliably producing both n-type and p-type 
Ti0 2 semiconducting materials are highly desirable. So far, doped Ti0 2 tends to yield only n- 
type semiconductors. However, it has recently been proposed that p-type Ti0 2 semiconducting 
materials may be obtained by nitrogen doping surface sites of low dimensional materials^. In this 
chapter we will discuss in detail the effects of quantum confinement and doping on the optical 
properties of Ti0 2 . 

Ti0 2 nanostructures are also one of the main components of hybrid solar cells. In a typical 
Gratzel cel^-, Ti0 2 nanoparticles with average diameters around 20 nm collect the photoelectron 
transferred from a dye molecule adsorbed on the surface 2 ^ 2 ^. Such processes are favoured by a 
proper energy level alignment between solid and organic materials, although the dynamic part of 
the process also plays an important role in the charge transfer. Clearly, Ti0 2 's characteristics of 
long quasi-particle lifetimes, high corrosion resistance, and relatively low cost, must be balanced 
with control of its energy level alignment with molecular states, and a fast electron injection at the 
interface. 

Despite all the engineering efforts, the main scientific goal remains to optimize the efficiency of 
solar energy conversion into readily available electricity. Different research approaches have been 
devoted to benefit from quantum size properties emerging at the nanoscal o 24 ' 25 , find an optimal 
donor- acceptor complex-2 4 - 2 ^, mix nanoparticles and one-dimensional structures, such as nanotubes 
or nanowires 2 ^ 2 !, and control the geometry of the Ti0 2 nano-assembly 2 ^ 

A clear theoretical understanding of Ti0 2 's optoelectronic properties is necessary to help un- 
ravel many fundamental questions concerning the experimental results. In particular, the properties 
of excitons, photo-injected electrons, and surface configuration in Ti0 2 nanomaterials may play a 
critical role in determining their overall behaviour in solar cells. For Ti0 2 at both the nanoscale 
and macroscale regime, it is necessary to first have a complete picture of the optical properties in 
order to clarify the contribution of excitons. 

Despite the clear importance of its surfaces and nanostructures, investigations of Ti0 2 bulk (see 
Q]) electronic and optical properties have not provided, so far, a comprehensive description of the 
material. Important characteristics, such as the electronic band gap, are still undetermined. Most of 
the experimental and computational work has been focused on synthesis and analysis of systems 
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FIG. 1: Unit cell (left), crystallographic structure (center) and Ti06 octahedrons (right) of rutile (top) and 
anatase (bottom). The lattice parameters in A are denoted a, b, and c, while R a x and R e q are the distances 
in A between a Ti 4+ ion and its nearest and next-nearest neighbour O 2 " ions, respectively. In the case of 
rutile the interstitial Ti impurity site is shown with a green circle (see left top). 

with reduced dimensionality. 

The experimental synthesis and characterization of nano structured materials is in general a 
costly and difficult task. However, predictions of a dye or nanostructure's properties from simula- 
tions can prove a great boon to experimentalists. Modern large-scale electronic structure calcula- 
tions have become important tools by providing realistic descriptions and predictions of structure 
and electronic properties for systems of technological interest. 

Although it will not be treated in this chapter, it is important to mention the problem of electron 
localization in reduced titania 2 ^. This will provide a glimpse of the complexity faced, from the 
theoretical point of view, when studying transition metal atoms. The localization of d electrons 
makes the accurate description of their exchange-correlation interaction a difficult task 30 . The 
electron localization in defective titania has been an open question from both experimental and 
theoretical points of view, and caused much controversy during the past few years^— . Oxygen 
vacancies are quite common in Ti0 2 , and their presence and behaviour can significantly affect 
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the properties of nanostructures. When an oxygen vacancy is created in Ti02, i.e. when Ti02 
is reduced, the two electrons coming from the removed O 2 " ion must be redistributed within the 
structure. One possibility is that these two extra electrons remain localized onto two Ti ions close 
to the O 2 " vacancy. In this way a pair of Ti 4+ ions become Ti 3+ ions. Another option is for the two 
extra electrons to delocalize along the whole structure, i.e. they do not localize on any particular 
Ti ion. Finally, an intermediate situation, with one electron localized and the other spread, is also 
possible. Concerning the Ti02 bulk, conventional density functional theory (DFT) calculations 
using either local density approximations (LDA) or generalized gradient approximations (GGA) 
for the exchange-correlation (xc)-functionals show a scenario with both electrons delocalized. On 
the other hand, hybrid functionals and Hartree-Fock calculations give rise to a situation with both 
electrons localized. For GGA+U calculations, the results are very sensitive to the value of the U 
parameter. For certain values of U both electrons remain localized, while for others there is an 
intermediate situation 3 ^^. 

Experimentally, there are electron paramagnetic resonance (EPR) measurements suggesting 
that the extra electrons are mainly localized on interstitial Ti 3+ ions 3 ^. These interstitial Ti 3+ ions 
are impurities placed at the natural interstices of the rutile lattice (seed)) and, similarly to the Ti 4+ 
ions of the pure lattice, they also form Ti0 6 octahedrons. Recent STM and PES experiments have 
shown that the interstitial Ti 3+ ions play a key role in the localization of the electrons when a bridge 
oxygen is removed from the Ti02 (110) surface 3 -^. These experiments concluded the controversial 
discussion about the localization of electrons in the bridge oxygen defective Ti0 2 (110) surface 
(see Refs.— ~— for more details). However, the problem remains unresolved for the bulk case. 

In summary, in this chapter we first analyze the full ab initio treatment of electronic and optical 
properties inland Dill before applying it to the two most stable bulk phases, rutile and anatase 
in HVl These are also the phases most easily found when nanostructures are synthesized. We will 
focus on their optical properties and excitonic behaviour. We then explore the possibility of tuning 
the oxide band gap using quantum confinement effects and dimensionality, by analyzing atomic 
clusters, nanowires and nanotubes in|V] A further component whose effect has to be evaluated is 
that of doping, which may further tune the optical behaviour by introducing electronic states in 
the gap, as presented in|VllCombining the effects of quantum confinement and doping is hoped to 
produce a refined properties control. In lVIIi we report some details on modeling for dye-sensitized 
solar cells, before providing a summary and our concluding remarks. 
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II. GROUND STATE PROPERTIES THROUGH DENSITY FUNCTIONAL THEORY 



DFT is a many-body approach, successfully used for many years to calculate the ground-state 
electronic properties of many-electron systems. However, DFT is by definition a ground state the- 
ory, and is not directly applicable to the study of excited states. To describe these types of physical 
phenomena it becomes necessary to include many-body effects not contained in DFT through 
Green's function theory. The use of many-body perturbation theory 3 ^, with DFT calculations as a 
zero order approximation, is an approach widely used to obtain quasi-particle excitation energies 
and dielectric response in an increasing number of systems, from bulk materials to surfaces and 
nanostructures. We present a short, general discussion of the theoretical framework, referring the 
reader to the books and the reviews available in the literature for a complete description (see, for 
instance, Refs. 39 and 40 ), before applying these methods to rutile and anatase Ti0 2 . 

As originally introduced by Hohenberg and Kohn (HK) in 1964 41 , DFT is based on the theorem 
that the ground state energy of a system of N interacting electrons in an external potential V ext {r) 
is a unique functional of the ground state electronic density. The Kohn and Sham 43 formulation 
demonstrates how the the HK Theorem may be used in practice, by self-consistently solving a set 
of one-electron equations (KS equations), 



where p(r) is the electronic charge density, (f>i(r) are the non-interacting KS wavefunctions, and 
K/7[p( r )] = ^ext( r ) + Vh[p(y)} + V xc [p(r)] is the effective one-electron potential. Here, Vh 
is the Hartree potential and V xc is the exchange-correlation potential defined in terms of the xc- 
functional E xc as V xc = J^y, which contains all the many-body effects. V xc is usually calculated in 
either LDA 4 ^ 4 - or GGA 4 ^ approximations. However, semi-empirical functionals are also available, 
called hybrids 4 -^, which somewhat correct the deficiencies of LDA and GGA for describing 
exchange and correlation. This is accomplished by including an exact exchange contribution. 

Other than the highest occupied molecular orbital (HOMO), KS DFT electronic levels do not 
correspond to the electronic energies of the many electron system. Indeed, the calculated KS band 
gaps of semiconductors and insulators severely underestimate the experimental ones. Experimen- 
tally, occupied states are accessible by direct photoemission, where an electron is extracted from 
the system (N — 1 ground state), while unoccupied states are accessible by inverse photoemission, 
where one electron is added to the system (N + 1 ground state). For isolated systems with a finite 
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(1) 
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number of electrons, the electronic gap may be obtained from the DFT calculated ground- state en- 
ergies with N + 1, N, and N — 1 electrons 4 ^. However, for periodic systems, adding or removing 
an electron is a non-trivial task. We therefore summarize in the following a rigorous method for 
describing excitations based on the Green's function approach. This method allows us to properly 
describe the electronic band structure. Further information and details about the Green's function 
approach may be found elsewher o 40 ' 49 ' — . 



III. EXCITED STATES WITHIN MANY BODY PERTURBATION THEORY 



When a bare particle, such as an electron or hole, enters an interacting system, it perturbs 
the particles in its vicinity. In essence, the particle is "dressed" by a polarization cloud of the 
surrounding particles, becoming a so-called quasi-particle. Using this concept, it is possible to 
describe the system through a set of quasi-particle equations by introducing a non-local, time- 
dependent, non-Hermitian operator called the self-energy S. This operator takes into account the 
interaction of the particle with the system via 

-~V 2 + V ext + V H * l (r,u) + J E(r,r',uj)^ i (r',u;)dr > = Eiiu^iir), (2) 

where ^(r) is the quasi-particle wavefunction. 

Since the operator £ is non-Hermitian, the energies Ei(u) are in general complex, and the 
imaginary part of £ is related to the lifetime of the excited particle 52 . The most often used ap- 
proximation to calculate the self-energy is the so-called GW method. It may be derived as the 
first-step iterative solution of the Hedin integral equations (see Refs.—,—, and^ 4 -), which link the 
Green's function G, the self-energy S, the screened Coulomb potential W, the polarization P and 
the vertex Y. In practice[2]is not usually solved directly, since the KS wavefunctions are often very 
similar to the GW ones 5 ^. For this reason it is often sufficient to calculate the quasi-particle (QP) 
corrections within first-order perturbation theory (GqWq) Moreover, the energy dependence 
of the self-energy is accounted for by expanding E in a Taylor series, so that the QP energies e i oWo 
are then given by 

£?0 W = £ KS + ZM KS n£ K S) _ V KS^KS^ (3) 

where Zi are the quasi-particle renomalization factors described in Refs.— and 4 ^. For a large num- 
ber of materials, the G W approximation of E works quite well at correcting the KS electronic 
gap from DFT. 
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Concerning optical properties, the physical quantity to be determined in order to obtain the opti- 
cal spectra is the macroscopic dielectric function €m(w). This may be calculated at different levels 
of accuracy within a theoretical ab initio approach. A major component in the interpretation of 
the optical measurements of reduced dimensional systems are the local-field effects (LFE). These 
effects are especially important for inhomogeneous systems. Here, even long wavelength exter- 
nal perturbations produce microscopic fluctuations of the electric field, which must be taken into 
account. However, LFE are also important for bulk phases such as anatase and rutile Ti02- They 
must be taken into account in the evaluation of optical absorption, and to calculate the screened 
interaction W used in GW. The effect becomes increasingly important when going to lower di- 
mensional systems. 

It is well known 55 that for inhomogeneous materials e M (cu) is not simply the average of the 
corresponding microscopic quantity, but is related to the inverse of the microscopic dielectric 
matrix by 

e M (a;) = lim— — -. (4) 

The microscopic dielectric function may be determined within the linear response theory^, the 
independent-particle picture by the random phase approximation (RPA), and using eigenvalues 
and eigenvectors of a one-particle scheme such as DFT or GW. There is also a different formu- 
lation which includes LFE in the macroscopic dielectric function. This becomes useful when the 
electron-hole interaction is included in the polarization function. This formulation allows us to 
include, via the Bethe-Salpeter equation (BSE), excitonic and LFE on the same footing. In so 
doing, inverting of the microscopic dielectric matrix is avoided. The complete derivation may be 
found in Appendix B of Ref. 57 . 

So far, in RPA, we have treated the quasi-particles as non-interacting. To take into account the 
electron-hole interaction, a higher order vertex correction needs to be included in the polarization. 
In other words, the BSE, which describes the electron-hole pair dynamics, needs to be solved. 
As explained in Ref. 57 , the BSE may be written as an eigenvalue problem involving the effective 
two-particle Hamiltonian 

u(ni,n 2 ),(n 3 ,n 4 ) = (T? _ ft )X X _ X( f _f ) 

iJ eic \ J -'ri2 J-'ni ) u ni,n3 u n2,ri4 \Jn2 J nil 

x J dr 1 dr' 1 dr 2 dr' 2 <j> ni (n) 0* n2 (r[) S(r 1; r' 1; r 2 , r' 2 ) 0; 3 (r 2 ) n4 (r' 2 ). 
The kernel H contains two contributions: v, which is the bare Coulomb interaction with the long 



range part corresponding to a vanishing wave vector not included and W, the attractive screened 
Coulomb electron-hole interaction. Using this formalism and considering only the resonant part 
of the excitonic Hamiltonian 57 , the macroscopic dielectric function may be expressed as 57 

2 



€ M (u) = 1 + limu(q; 



E 



(v,c;k) 



q^O ^ (E X - W) 



(5) 



In |5] the dielectric function, differently from the RPA approximation, is given by a mixing of 
single particle transitions weighted by the excitonic eigenstates A\. These are obtained by the 
diagonalization of the excitonic Hamiltonian. Moreover, the excitation energies in the denominator 
are changed from e c — e v to Ex- The electronic levels are mixed to produce optical transitions, 
which are no longer between pairs of independent particles. The excitonic calculation is in general, 
from the computational point of view, very demanding because the matrix to be diagonalized may 
be very large. The relevant parameters which determine its size are the number of k-points in the 
Brillouin zone, and the number of valence and conduction bands, N v and N c respectively, which 
form the basis set of pairs of states. 

Calculations performed for insulators and semiconductors show that the inclusion of the 
electron-hole Coulomb interaction yields a near-quantitative agreement with experiment. This 
is not only true below the electronic gaps, where bound excitons are generally formed, but also 
above the continuum edge. The same results apply to the titania-based materials investigated here, 
as shown in the following section. 



IV. THE BULK PHASES OF TI0 2 : ROLE OF MANY BODY EFFECTS 

Despite the importance of its surfaces and nanostructures, the most recent measurements of 
TiCVs bulk (see Fig. 1) electronic and optical properties were performed in the 1960s, with a few 
exceptions. Here we aim to review the existing results obtained using a variety of experimental 
techniques and ab initio calculations, in order to elucidate the known properties of Ti0 2 . Previous 
data will be compared with a complete, consistent ab initio description, which includes many body 
effects when describing electronic and optical properties. Most experimental and theoretical data 
reported refers to the rutile phase, while anatase in general has been less studied. However, anatase 
has received more attention recently, due in part to its greater stability at the nanoscale compared 
to rutile. 
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As we will see, while a general agreement seems to exist concerning the optical absorption 
edge of these materials, values for basic electronic properties such as the band gap still have a 
large degree of uncertainty. 

The electronic properties of valence states of rutile Ti0 2 have been investigated experimentally 
by angle-resolved photoemission spectroscopy 58 , along the two high symmetry directions (A and 
S) in the bulk Brillouin zone. The valence band of Ti0 2 consists mainly of O 2p states partially 
hybridized with Ti 3d states. The metal 3c? states constitute the conduction band, with a small 
amount of mixing with O 2p states. This photoemission data was compared to calculations per- 
formed with both pseudopotentials and linear muffin-tin orbital (LMTO) methods, which gave a 
direct gap of 2 eV in both cases. On the other hand, within the linear combination of atomic 
orbitals (LCAO) method, a gap of 3 eV was obtained. 

From the symmetry of the TiC>6 octahedrons (see [D, d states are usually grouped into low 
energy t 2g and high energy e g sub-bands. It is important to note that, from ultraviolet photoe- 
mission spectroscopy (UPS) data, it has been deduced that the electronic gap for rutile is at least 
4 eV. This is the observed binding energy of the first states below the Fermi energy. This is in 
agreement with previous reported data from electron energy loss spectroscopy 5 ^, and from other 
UP$ffi~— results. The electronic structure of rutile bulk has also been described using other exper- 
imental techniques, such as electrical resistivity 63 , electroabsorption 6 ^ 5 -, photoconductivity and 
photoluminescenc o 66 ' 67 , X-ray absorption spectroscopy (XAS)^— , resonant Raman spectra 6 - 7 -^ 
photoelectrochemical analysi s 75 ' 76 and UPS 77 . All of these experiments have provided many im- 
portant details of its electronic properties, in particular concerning the hybridization between Ti 
3d and O 2p states. However, the electronic band gap, corresponding to the difference between the 
valence band maximum (VBM) and the conduction band minimum (CBM), has not been obtained 
directly from any experimental data. Although the electronic band gap could be measured using 
combined photoemission and inverse photoemission experiments, such experiments do not appear 
in the literature. 

The same discussion is valid for the anatase crystalline phase. Even though there are several 
XAS measurements concerning its electronic structur o 73i78]79 , photoemission data is completely 
lacking for anatase. In the absence of more recent and refined experimental results for rutile, and 
due to the lack of results for anatase, we are left with an estimate of 4 eV for the electronic gap of 
rutile TiC>2. It is this value which we shall use as a reference in the following discussion. 

We will now review the experimental results for optical properties of Ti0 2 . Such measurements 
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are of great interest for the photocatalytic and photovoltaic applications of this material. From the 
optical absorption spectra of both phases 80 , the room temperature optical band gap is found to be 
3.0 eV for rutile, and 3.2 eV for anatase. 

The absorption edge has been investigated in detail for rutile by combining absorption, photo- 
luminescence, and Raman scattering techniques 81 . These techniques provide a value for the edge 
of 3.031 eV associated to a 2p xy exciton, while a lower energy Is quadrupolar exciton has been 
identified below 3 eV. The first dipole allowed gap is at 4.2 eV 81 according to the combined results 
of these three techniques. Concerning anatas o 82 ' 83 , the optical spectrum have been recently re- 
evaluated 8 -, confirming the 3.2 eV value for the edge. The fine details of anatase's spectrum have 
also been recently investigated 84 . Data on the Urbach tail has revealed that excitons in anatase are 
self-trapped in the octahedron of coordination of the titanium atom. This is in contrast to the rutile 
phase, where excitons are known to be free due to the different packing of rutile's octahedra.— 

In general, measurements of optical properties can be significantly affected by the presence of 
defects, such as oxygen vacancies, and by phonons. Both defects and phonons will be present 
in any experimental sample of the material at finite temperatures. These observations have to 
be kept in mind when directly comparing experimental measurements with the theoretical results 
presented in the following. Moreover, there is a general trend in theoretical-computational studies 
to compare the theoretical electronic band gaps with the experimental optical gap value o 29i85i85 ~ — 
derived from the above mentioned experiments. 

It should be remembered that almost by definition, the optical gap is always smaller than the 
electronic band gap. This is because the two types of experiments (photoemission, and optical 
absorption) provide information on two different physical quantities. Reverse photoemission ex- 
periments involve a change in the total number of electrons in the material (N — > N + 1), while 
optical absorption experiments do not (N — > N*). The latter involves the creation of an electron- 
hole pair in the material, with the hole stabilizing the excited electron. For this reason, comparison 
between experimental and theoretical data, and the resulting discussion, must take into account 
the proper quantities. 

The theoretical investigations presented in the literature of TiC^'s structural, electronic, and 
optical properties are at varying levels of theory and thus somewhat inconsistent. A comprehensive 
description of properties of both phases in the same theoretical and computational framework is 
still missing. Here we present, in a unified description and by treating with the same method for 
the two phases, the electronic and optical properties of the two most stable phases of bulk titania. 
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The ab intio calculations performed yield results 9 in quite good agreement with the few available 
pieces of experimental data. 

A combination of DFT— and many body perturbation theory (MBPT) methods is a reliable 
and well established toolkit to obtain a complete analysis of electronic and optical properties for 
a large class of materials and structures. In this DFT + GW + BSE framework, the properties of 
the two bulk phases of TiC>2 may be properly analyzed. Their structural and energetic properties 
have been calculated 89 using DFT, as a well established tool for the description of ground state 
properties. However, the DFT gap is, as expected, significantly smaller than the experimental gap, 
with the relative positions of the s, p, and d levels also affected by this description. To address this, 
standard G W calculations may be applied to obtain the quasi-particle corrections to the energy 
levels, starting from DFT eigenvalues and eigenfunctions. Finally the electron-hole interaction is 
included, to properly describe the optical response of the system. 

The description of ground state propertie s 86i9Q performed in the framework of DFT are gen- 
erally quite good, with the structural description of TiC>2 systems in reasonable agreement with 
experiments 89 . The lattice constants are within 2% of experiment, while bulk modulii are within 
10% of the experimental results 86 ' 90 , as is often found for DFT. However, DFT incorrectly pre- 
dicts the anatase phase to be more stable than the rutile one, even for a small energy difference, 
independently of the xc-functional used 86 . 

Even if DFT is not an excited state method, the KS wavefunctions are often used to evaluate the 
band structure along the high symmetry directions (cf. [2] for rutile), the density of states, and the 
spatial behaviour of wavefunctions involved in the relevant bonds in the system 5818919 ' 192 . The KS 
electronic gap, corresponding to the difference between the VBM and the CBM, is 1 .93 eV and 
2.16 eV for rutile (direct gap) and anatase (indirect gap), respectively. These are underestimations 
by almost 2 eV of the available experimental data 5 ^,. However, the overall behaviour of the band 
dispersion of KS levels is reasonable, with valence bands mainly given by O 2p states, and Ti 3d 
states forming the conduction bands. 

The application of standard GW methods gives gaps of 3.59 and 3.97 eV for rutile and 
anatase^, respectively. The value for rutile is again smaller than the one given by the UPS es- 
timation 5 ^, but still close to the experimental value of 4 eV. 

There exist a number of theoretical works, with calculations performed at different levels of 
DFT or including MBPT descriptions, for the electronic gap of rutile and anatase Ti0 2 85l85l87l91 ~ — . 
Therefore in literature it is possible to find for the electronic gap a quite large range of possible 
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FIG. 2: Electronic band structure of rutile bulk, along the high symmetry directions of the irreducible 
Brillouin zone, from GGA calculations ( ). and including the GoWq correction ( ). 

values, attributed to the gap of titanium dioxide, which are often erroneously compared with the 
experimental optical gap. 

The DFT-GGA values calculated^ are comparable to the ones obtained with a variety of dif- 
ferent DFT approaches, with different functionals, by using plane waves or localized basis meth- 
ods, and all-electron or pseudopotential approaches. Only the hybrid PBEO 46 and B3LYP^ xc- 
functionals give larger values. 

From quasi-particle calculations, the electronic gap of anatase has been estimated to be 
3.79 eV— by G W . However the more refined computational approach, because of its inclusion 
of a self-consistent evaluation of GW, yields a gap of 3.78 eV for rutile TiC^.^^ 

Moving to optical properties, and by applying the RPA method to both KS and QP energies, 
we obtain spectra (Fig. |3]) that do not in overall behaviour agree with experiment. Differences are 
clear both in absorption edge determination, and in the overall shape of the spectra. The inclu- 
sion of quasi-particle corrections at the GW level yields a rigid shift of the absorption spectrum, 
moving the edge at higher energies, due to the opening of the gap. However, the shape of the 
absorption is quite unaffected, because the interaction is still treated using an independent quasi- 
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FIG. 3: Imaginary part of the dielectric constant for rutile (left), and anatase (right), in-plane polarization, 

calculated by GGA RPA ( ), using Go Wo on top of GGA ( ), and via the Bethe-Salpeter equation 

(BSE) ( ). The experimental spectrum ( — • — ) from Refs.— and^ 2 - is also shown for comparison. 

particle approximation. A substantially better agreement may be obtained^^ by solving the BSE, 
which takes into account both many body interactions and excitonic effects. Indeed, it produces a 
good description of absorption spectra and excitons, as shown in Fig. [3] 

The optical absorption spectra calculated for the two phases, with polarization along the x- 
direction of the unit cells, are provided in [3] The spectra given by independent-particle transitions 
present two characteristic features. First, the band edge is underestimated, due to the electronic 
gap underestimation in DFT. Second, the overall shape of the spectrum is, for both phases, and 
both orientations, different from the experiment, in the sense that the oscillator strengths are not 
correct. The inclusion of the quasi-particle description, which should improve the electronic gap 
description, does not improve the overall shape of the spectrum. The absorption edge is, however, 
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FIG. 4: Spatial distribution (yellow isosurfaces in arbitrary units) of the partial dark, dark, and optically 
active first three excitons in rutile. The hole position is denoted by the light green dot. 

shifted at higher energies, even higher than expected from experiments. The description of optical 
properties within the two interacting quasi-particle approximation (by solving the BSE) definitely 
improves the results. The absorption edge is now comparable to the experimental one, with the 
optical gap estimated from our calculations in good agreement with the available data. Further, the 
shape of the spectrum is now well described, with a redistribution of transitions at lower energies. 
The agreement is generally good for both phases 89 . 

The nature of the exciton is still under debate in Ti0 2 materials^. The experimental binding 
energy is of 4 meV and some uncertainty exists for the exact determination of optical edge. More- 
over, the exciton is localized in one of the two phases, and delocalized in the other one, at least 
based on experimental results. However, an explanation for this behaviour is so far missing. Re- 
fined measurements^ give an exciton of 2p xy character at 3.031 eV, and a Is quadrupolar exciton 
at energy lower than 3 eV. From ab initio calculations with x polarization, two dark (optically 
inactive) or quasi dark excitons are located in rutile at 3.40 and 3.55 eV (denoted by E B in 
IU). At the same time, the optically active exciton, with 4 meV of binding energy, is located at 
3.59 eV (Ec in|4]). The spatial distribution of the first three excitons is plotted in|4] The transitions 
are from O 2p states to Ti 3d states of the triplet t 2 g> as expected. While the first two optically 
forbidden transitions involve Ti atoms farther away from the excited O atom, the optical active 
transition involves states of the nearest neighbour Ti atoms. 

In this section we have endeavoured to clarify though the application of a consistent description, 
the properties of the two main crystalline phases of titania. Particular attention has been taken to 
how the inclusion of a proper description of exchange and correlation effects can improve the 
description of both electronic and optical properties of T1O2. 
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Now that bulk properties are known, from a theoretical point of view, at the level that the 
state-of-the-art ab initio techniques allow us to reach, we can attempt to describe how quantum 
confinement induced by reduced dimensionality, and doping both effect the properties of Ti0 2 . 
Our final aim is to demonstrate how we may tune the electronic and optical gap of nanostructures 
for photocatalytic and photovoltaic applications. 

V. USING NANOSTRUCTURE TO TUNE THE ENERGY GAP 

Having described in the previous section the electronic structure of the two bulk Ti0 2 phases, 
we will now turn our attention to the influence of nanostructure on the energy gap. This is an 
area which has received considerable attention in recent years^— , in part due to the inherently 
high surface to volume ratio of nanostructures. It is hoped that this will allow materials with 
shorter quasi-particle lifetimes to function effectively for photocatalytic activities, since excitons 
are essentially formed at the material's surface. However, this advantage is partly countered by 
quantum confinement effects, which tend to increase the energy gap in nanostructures. These 
competing factors make the accurate theoretical determination of energy gaps in nanostructured 
materials a thing of great interest. 

However, as shown in the previous section, standard DFT calculations tend to underestimate 
electronic band gaps for bulk Ti0 2 by approximately 2 eV, due in part to self-interaction error s 96 i 91 . 
These errors arise from an incomplete cancellation of the electron's Coulomb potential in the 
exchange-correlation (xc)-functional. 

This may be partially addressed by the use of hybrid functionals such as B3LYP— , which 
generally seem to improve band gaps for bulk system s 32187 ' 98 . However, such calculations are com- 
putationally more expensive, due to the added dependence of the xc-functional on the electron's 
wavefunction. Moreover, B3LYP calculations for Ti0 2 clusters largely underestimate the gap rel- 
ative to the more reliable difference between standard DFT calculated ionization potential I p and 
electron affinity E a . For isolated systems such as clusters, the needed energetics of charged species 
are quantitatively described by standard DFT. Also, B3LYP and RPBE" calculations provide the 
same qualitative description of the trends in the energy gaps for Ti0 2 , as seen in [5] Another 
methodology is thus needed to describe the band gaps of periodic systems. 

GqWq is probably the most successful and generally applicable method for calculating quasi- 
particle gaps. For clusters it agrees well with I p — E a , and it has been shows to produce reliable 
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FIG. 5: Energy gap e gap in eV versus Ti02 structure for OD (Ti02) n clusters (n < 9), ID Ti02 (2,2) 
nanorods, (3,3) nanotubes, (4,4) nanotubes, 2D HexABC and anatase layers, and 3D anatase surface, anatase 
bulk, and rutile bulk phases. DFT calculations using the highest occupied and lowest unoccupied state 
gaps with the standard GGA RPBE xc-functional (O) and the hybrid B3LYP xc-functional (T), are com- 
pared with DFT I p — E a calculations (•), GqWq quasi-particle calculations (►), and experimental results 
(♦) 6 1 11 i 13 1 20 i 94 . Schematics of representative structures for each dimensionality are shown above and taken 
from Ref . 1 1 . 

results for bulk phases. On the other hand, G W calculations describe an N — > N + 1 transition 
where the number of charges is not conserved, rather than the electron-hole pair induced by pho- 
toabsorption, which is a neutral process. Indeed, a description in terms of electronic gap cannot 
be compared with, or provide direct information on the optical gap, which is the most investigated 
quantity, due to its critical importance for photocatalytic processes. 

[5] shows that for both 3D and 2D systems, RPBE gaps underestimate the experimental optical 
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gaps by approximately 1 eV. For ID and OD systems, there is a much larger difference of about 
4 eV and 5 eV respectively, between the RPBE gaps and the I p — E a and GqWq results. This 
increasing disparity may be attributed to the greater quantum confinement and charge localization 
in the ID and OD systems, which yield greater self-interaction effects. The B3LYP gaps also tend 
to underestimate this effect, simply increasing the RPBE energy gaps for both OD and 3D systems 
by about 1.4 eV. 

On the other hand, the RPBE gaps reproduce qualitatively the structural dependence of the 
I p — E a , GoWo, and experimental results for a given dimensionality, up to a constant shift. This is 
true even for 3D bulk systems, where standard DFT does not correctly predict rutile to be the most 
stable structure^. 

To summarize, quantum confinement effects seem to increase the energy gap significantly for 
both OD and ID systems, while 2D and 3D systems may be much less affected. This suggests that 
2D laminar structures are viable candidates for reducing the minimal quasi-particle, while leaving 
the band gap nearly unchanged. However, a more accurate description of the photoabsorption 
properties of these novel nanostructures, perhaps using the methodologies recently applied to bulk 
TiOy^, still remains to be found. 

VI. INFLUENCE OF BORON AND NITROGEN DOPING ON TI0 2 'S ENERGY GAP 

The doping of Ti0 2 nanostructures has received much recent attention as a possible means 
for effectively tuning Ti0 2 's band gap into the visible range^^— . Recent experiments suggest 
substituting oxygen by boron or nitrogen in the bulk introduces mid-gap states, allowing lower 
energy excitations. However, to model such systems effectively requires large supercells, both 
to properly describe the experimental doping fractions of < 10%, and to ensure dopant-dopant 
interactions are minimized. 

[6] shows the DFT calculated DOS and structures for the most stable boron doped and nitrogen 
doped Ti0 2 (2,2) nanorods. The highest occupied state is also shown as isosurfaces of ±0.05e/A 3 
in the side views of the doped structures. 

As with Ti0 2 clusters, the influence of boron dopants on Ti0 2 nanorods may be understood in 
terms of boron's weak electronegativity, especially when compared with the strongly electronega- 
tive oxygen. Boron prefers to occupy oxygen sites which are 2-fold coordinated to neighbouring 
titanium atoms. However, as with the OD clusters, boron's relatively electropositive character 
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FIG. 6: Total density of states (DOS) in eV" 1 per T1O2 functional unit vs. energy e in eV for undoped ( ), 

boron doped ( ) and nitrogen doped ( ) (2,2) T1O2 nanorods from standard DFT GGA RPBE xc- 

functional and (inset) GqWq quasi -particle calculations. The highest occupied states for boron and nitrogen 
doped (2,2) Ti02 nanorods are depicted by isosurfaces of ±0.05e/A 3 in the structure diagrams to the left. 

induces significant structural changes in the ID structures, creating a stronger third bond to a 
neighbouring three-fold coordinated oxygen via an oxygen dislocation, as shown in[6] This yields 
three occupied mid-gap states localized on the boron dopant, which overlap both the valence band 
O 2p n and conduction band Ti 3d xy states, as shown in [6l Boron dopants thus yield donor states 
near the conduction band, which may be photocatalytically active in the visible region. However, 
the quantum confinement inherent in these ID structures may stretch these gaps, as found for the 
GqWo calculated DOS shown in the inset of [6l 

On the other hand, nitrogen dopants prefer to occupy oxygen sites which are 3-fold coordinated 
to Ti, as was previously found for the rutile Ti0 2 surface^^. This yields one occupied state at the 
top of the valence band and one unoccupied mid-gap state in the same spin channel. Both states 
are localized on the nitrogen dopant but overlap the valence band O 2p n states, as shown in [6] 
Nitrogen dopants thus act as acceptors, providing localized states well above the valence band, as 
is also found for the GqWq calculated DOS shown in the inset of [6l 

Although nitrogen dopants act as acceptors in Ti0 2 ID structures, such large gaps between 
the valence band and the unoccupied mid-gap states would not yield p-type semiconductors. This 
may be attributed to the substantial quantum confinement in these ID structures. However, for 
2D and 3D systems, it is possible to produce both p-type and n-type classical semiconductors, as 
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discussed in Ref. 11 . This has recently been shown experimentally in Ref. 19 , where co-doping of 
anatase Ti02 with nitrogen and chromium was found to improve the localization of the acceptor 
states, and reduce the effective optical gap. By replacing both Ti 4+ and O 2 " atoms with dopants 
in the same Ti0 6 octahedral, it should be possible to "tune" the optical band gap to a much finer 
degree. 

Whether calculated using RPBE, I p — E a or GqWq, the energy gaps for both boron and nitrogen 
doped Ti0 2 nanostructures are generally narrowed, as shown inf^a) and (b). However, for nitrogen 
doped (Ti02) n clusters where nitrogen acts as an acceptor (n = 5, 6, 9), the energy gap is actually 
increased when spin is conserved, compared to the undoped clusters in RPBE. This effect is not 
properly described by the N — > N + 1 transitions of I p — E a , for which spin is not conserved for 
these nitrogen doped clusters. On the other hand, when nitrogen acts as a donor (n = 7, 8) the 
smallest gap between energy levels does conserve spin. 

The boron doped Ti0 2 nanorods and nanotubes have perhaps the most promising energy gap 
results of the Ti0 2 structures, as seen in [3a)- Boron dopants introduce in the nanorods localized 
occupied states near the conduction band edge in both RPBE (cf. [6]) and GoWq (cf. inset of [6]) 
calculations. On the other hand, nitrogen doping of nanorods introduces well defined mid-gap 
states, as shown in [6] However, to perform water dissociation, the energy of the excited electron 
must be above that for hydrogen evolution, with respect to the vacuum level. This is not the case 
for such a mid-gap state. This opens the possibility of a second excitation from the mid-gap state 
to the conduction band. However, the cross section for such an excitation may be rather low. 

For boron doping of 2D and 3D structures, the highest occupied state donates its electron almost 
entirely to the conduction band, yielding an n-type semiconductor. Thus at very low temperatures, 
the RPBE band gap is very small. The same is true for n-type nitrogen doped bulk anatase. 
For these reasons the energy gap between the highest fully occupied state and the conduction 
band, which may be more relevant for photoabsorption, is also shown. These RPBE gaps are still 
generally smaller than those for their undoped Ti02 counterparts, shown in [51 

In summary, for both boron and nitrogen doped clusters we find RPBE gaps differ from I p — E a 
by about 3 eV, while for nitrogen doped anatase the RPBE gap differs from experiment by about 
0.6 eV. Given the common shift of 1 eV for undoped 2D and 3D structures, this suggests that both 
boron and nitrogen doped 2D Ti0 2 structures are promising candidates for photocatalysis. Further, 
the boron and nitrogen doped ID nanotube results also warrant further experimental investigation. 
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FIG. 7: Influence of doping on the energy gap e gap in eV versus T1O2 structure for OD (Ti02) n clusters 
(n < 9), ID T1O2 (2,2) nanorods, (3,3) nanotubes, (4,4) nanotubes, 2D HexABC and anatase layers, and 
3D anatase surface and bulk. The energy gaps for (a) boron doped systems from DFT calculations using 
the highest occupied and lowest unoccupied state gaps with the standard GGA RPBE xc-functional (A) are 
compared with DFT I p — E a calculations (A), and GqWq quasi-particle calculations (►), and (b) nitrogen 
doped systems from DFT calculations using the highest occupied and lowest unoccupied state gaps with the 
standard GGA RPBE xc-functional (□) are compared with DFT I p — E a calculations (■), and Go Wo quasi- 
particle calculations (►), and experimental (♦) results ^ 1 M 3 ! 20 ? 94 , Small open symbols denote transitions 
between highest fully occupied states and the conduction band. 

VII. SOLAR CELLS FROM TI0 2 NANOSTRUCTURES : DYE-SENSITIZED SOLAR CELLS 



TiC>2 is so far the most widely used solid material in the development of solar cell devices 
based on hybrid architectures^^. In these devices, the dye, synthetic or organic, absorbs light, 
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FIG. 8: Absorption spectrum calculated by TD Adiabatic LDA of the Ru-dye N719, whose structure is 
shown in the inset. 

and electrons excited by the phonons are injected into the underlying oxide nanostructure. The 
hybrid system must therefore satisfy several requirements: (1) a proper absorption range for the 
dye, (2) a fast charge transfer in the oxide, (3) a slow back-transfer process, and (4) an easy 
collection and conduction of electrons in the oxide. 

While absorption properties may be easily tuned at the chemical level by changing or adding 
functional groups, a critical point is to understand, and therefore control, the process at the in- 
terface. Indeed, simply having a good energy level alignment is not sufficient because the fast 
electron injection process is dynamical. Since the experimental characterization of complex sys- 
tems (networks of nanostructures, with adsorbed dyes, and in solution) is quite complicated, the 
theoretical description of such systems can be of fundamental importance in unravelling the pro- 
cesses governing the behaviour of dye-sensitized solar cells (DSSC). 

The most popular technique for studying these dynamical processes is time dependent DFT 
(TDDFT). For a detailed discussion of the success and possible limitations of this method when 
applied to hybrid systems, we refer the reader to Ref.— . TDDFT is a generalization of DFT which 
allows us to directly describe excited states. For example, it has previously been used successfully 
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to calculate the optical absorption of large organic molecules ®, such as Ru-dyes, or indolines. 
This same technique has also been applied recently to hybrid systems used for photovoltaic ap- 
plications. We want here to highlight that the charge injection transfer has also been modelled 
for molecules adsorbed on Ti0 2 clusters^—, giving estimations as fast as 8 fs for the injection 
time 107 . TDDFT is therefore a powerful tool to understand dynamical electronic effects in systems 
as large and complex as hybrid organic-oxide solar cells. 

VIII. CONCLUSION 

In this chapter we showed how an oxide of predominant importance in nanotechnological and 
environmental fields, such as photocatalysis and photovoltaic s, can be successfully investigates 
by state-of-the-art ab initio techniques. The optical and excitonic properties of the bulk phases 
can be properly described by MBPT techniques. We also showed how the electronic properties 
of TiC>2 may be "tailored" using nano structural changes in combination with boron and nitrogen 
doping. While boron doping tends to produce smaller band gap n-type semiconductors, nitrogen 
doping produces p-type or n-type semiconductors depending on whether or not nearby oxygen 
atoms occupy surface sites. This suggests that a p-type TiC>2 semiconductor may be produced 
using nitrogen doping in conjunction with surface confinement at the nanoscale. We also showed 
that it has been proved how, in the field of photovoltaics, TDDFT is a powerful tool to understand 
the mechanism of charge injection at organic-inorganic interfaces. 
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